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The properties of surfaces with charge-regulated patches are studied using non-linear Poisson- 
Boltzmann theory. Using a mode expansion to solve the non-linear problem efficiently, we reveal 
the charging behaviour of Debye-length sized patches. We find that patches charge up to higher 
charge densities if their size is relatively small and if the patches are well separated. The numerical 
results are used to construct a basic analytical model which predicts the average surface charge 
density on surfaces with patchy chargeable groups. 



I. INTRODUCTION 

Most surfaces that arc immersed in an aqueous 
solution obtain a net charge due to ion adsorption or 
dissociative processes at the surface. The resulting 
electrostatic force between such surfaces is essential 
for understanding the stability, osmotic pressure and 
ffocculation behaviour of colloidal suspensions. In 
the standard (linear) screening picture, like-charged 
surfaces repel each other at distances of the order of the 
Debye screening length, [TH3] due to overlapping clouds 
of screening ions in the vicinity of the two surfaces. 
This electrostatic repulsion, combined with short-range 
attractive Van der Waals forces if there is a dielectric 
contrast between the colloidal particles and the solvent, 
is a basic result of Derjaguin-Landau-Verwey-Overbeek 
(DLVO) theory which dates back to the 1940's.[TH3| More 
recent studies, based on nonlinear Poisson-Boltzmann 
theory, also find strictly repulsive electrostatic forces 
between pairs of like-charges surfaces. [1HH| Nevertheless, 
there are also experimental reports of attraction between 
like-charge colloids at ranges much longer than those 
of Van der Waals forces. Ion-ion correlations, which 
are ignored in the mean-field type Poisson-Boltzmann 
theory, might be an explanation for these observed 
attractions in the case of multi-valent ions. |10IfT2] 
However, evidence for electrostatic attractions has also 
been reported for suspensions with only monovalent 
ions. |13lH5] causing heated debates in the literature on 
the break-down of the classic DLVO theory due to many- 
body effects, the vicinity of glass walls, hydrodynamic 
forces, etc. Interestingly, it has also been suggested 
that charge inhomogeneities can be responsible for 
these attractions, P2H20] where the heterogeneity of the 
surface charge may be due to an incidentally present or 
purposely designed underlying chemical structure, or by 
clustering of adsorbed surfactants. 

Apart from the ill-understood electrostatic attractions 
in some systems, another good reason for considering het- 
erogeneously charged surfaces in more detail stems from 
the fascinating recent advances in the chemical synthe- 
sis of a large class of novel patchy nanoparticles, fea- 
turing not only corners, edges, and facets due to their 
finite size, but also spots or stripes. 21 Understanding 



the large-scale self-assembly properties of these newly 
available nanoparticles is an important ongoing scientific 
quest that requires effective particle-particle interactions 
as an input. For this reason a better understanding is 
needed of the relations between the chemical heterogene- 
ity of patchy particles, the resulting surface charge den- 
sity, and the nature and geometry of the ionic screening 
cloud that ultimately dictates the effective electrostatic 
interactions. In this article we explore some of these re- 
lations within nonlinear Poisson-Boltzmann (PB) theory 
in the relatively simple geometry of chargeable stripes on 
a planar surface in contact with a bulk electrolyte. 

In PB theory a key role is being played by the bound- 
ary conditions (BC's), in particular those on the surface 
between the suspending electrolyte and the suspended 
colloidal nanoparticles. The most common BC is to pre- 
describe the surface charge density of the colloidal par- 
ticle, thereby imposing a fixed discontinuity of the dis- 
placement field at both sides of the surface while the sur- 
face potential itself increases upon the approach of an- 
other like-charged surface. This type of constant-charge 
BC was shown to be realistic for e.g. strongly acidic 
homogeneous surfaces that ionise completely in a po- 
lar solvent such as water. [22] An alternative is to pre- 
describe the electrostatic potential on the surface, such 
that the particles can adjust their charge density if an- 
other surface approaches. This constant-potential BC 
has turned out to be realistic for surfaces not too far 
from a point of zero charge. [53] It is not clear, how- 
ever, which of these BC's is realistic for surfaces with 
a heterogeneous chemical composition. For instance, it 
is not trivial how the charge of a highly charged patch 
on an otherwise weakly charged or neutral surface is 
distributed, and how the induced electrostatic potential 
propagates to affect the charging of nearby surface groups 
of other patches. The problem requires us to treat the 
interplay between the electrostatic potential and the sur- 
face charge density at the more microscopic level of a 
chemical equilibrium of attaching and detaching ions as 
modeled by charge-regulation BC's. [23 25 Such an ap- 
proach has been applied to fit force measurements be- 
tween two homogeneously charged plates. [2"3T - [2"5] or to 
compute forces between periodically modulated charged 
plates within linearised PB theory. [26] Here, however, we 
combine charge regulation BC's with nonlinear PB the- 
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L = A + D 2 



in dielectric medium gives rise to a boundary condition 




Figure 1: Pictorial representation of the model. The striped 
patches are shaded to distinguish these from the rest of the 
plate and will be referred to as region "1". The area in be- 
tween the patches will be referred to as region "2" . 



ory for a single spatially heterogenous surface. Using a 
newly developed numerical scheme we will reveal how 
surfaces with finite-sized discrete patches |2"6"H2"5] tend to 
charge up and get screened, where the full non-linear PB 
theory allows us to deal with highly charged patches next 
to oppositely charged or neutral areas with strongly vary- 
ing ion concentrations perpendicular and parallel to the 
surface. 



II. THEORY 

We consider a solid medium in the half space z < 0, 
characterized by a dielectric constant ep. The surface of 
this medium at z = is considered to be chemically het- 
erogeneous, e.g. with a stripe pattern such as depicted in 
Fig.[l] The half space z > is a bulk solvent with dielec- 
tric constant e and volume V at temperature T. This sol- 
vent is assumed to be in thermal and diffusive equilibrium 
with a reservoir at z — > oo that contains point-like mono- 
valent cations and anions, both at a concentration p s . 
Since treating the ions in a mean-field fashion is in most 
of the cases an accurate approach in the case of monova- 
lent ions in water or the more polar (higher e) oils. [291 [50] 
we describe the distributions of ions for z > by the 
Boltzmann distributions p±(r) = p s exp(=F$(r)). Here 
/cBT$(r)/e is the electrostatic potential at r = (x,y,z), 
with e the elementary charge and fceT the Boltzmann 
constant. We assume that §(x,y,z —> oo) = in the 
reservoir, and that the densities of both ion species are 
zero for z < due to hard-core repulsions. The Poisson 
equation relates the charge density to the Laplacian of 
the potential, V 2 <I>(r) = — 47tAb (p+ (r) — p_(r)), where 
we introduce the Bjerrum length Ar, = e 2 /ek^T. What 
follows is the Poisson-Boltzmann (PB) equation, 



V 2 $(r) = 



k 2 sinh$(r) z > 0; 







z < 0, 



(1) 



where the Debye screening length is given by = 
(87tAbPs) -1 / 2 . At the fluid-solid interface, the presence 
of a surface charge density ea(x, y) and a stepwise change 



V$(x,y,0 + ) ■ z = -4nX B cr(x,y) 

+ —V$(x,y,0~) ■ z. 



(2) 



Other boundary conditions ensure that the potential is 
continuous at z = 0, vanishes at infinity, and require that 
the electric field in the medium vanishes far away from 
the interface: 

$(z,y,0+) = $(x,y,0-); (3) 
lim $(x,y,z) = 0; (4) 

z— >oo 

lim $'(x,y,z) = 0. (5) 

z— > — oo 

Here, and below, a prime denotes a partial derivative 
w.r.t. the z-coordinate. 



A. Charge regulation 

The majority of models which apply PB theory assume 
either a fixed charge density or a fixed potential at sur- 
faces. In reality the charging is often regulated; neither 
the charge density nor the surface potential is fixed, since 
both depend on the local density of ions in the surround- 
ing liquid. There is a subtle interplay as the local surface 
charge density depends on the local surface potential and 
vice versa. 

We consider 4 major charge generation mechanisms, 
involving cationic/anionic adsorption/desorption, 



S + H+ *=♦ SH+ with K$ = 



SB T± S+ + B" with 



S + B" 



SH 



SB" with K A = 



H+ with K D = 



[S][H+] 
[SH+] ' 
[S + ][B-] 
[SB] ' 
[S][B-] 
[SB"] ' 
d [S~][H + ] 



[SH] ' 



(6a) 
(6b) 
(6c) 
(6d) 



where S denotes a surface group, H + represent a cation, 
and B~ an anion, with concentrations [S], [H + ], and [B~] 
in the vicinity of the surface. Eqs. (6a) and ([6c]) describe 



associative charging in which an ion from the liquid is ad- 
sorbed at the surface. Eqs. (6b I and (6d) describe disso- 



ciative charging, where ions are released from the surface 
into the liquid. These chemical reactions are character- 
ized by reaction constants K± D . We assume the charg- 
ing to occur at discrete sites at the surface, and define 
a surface site density m±(x,y) and m±(x,y) for asso- 
ciatively and dissociatively chargeable sites respectively. 
Since each surface site is either charged or neutral, we 
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have 



[S] + [SH+] 
[S+] + [SB] 
[S] + [SB-] 
IS"] + [SH] 



m A (x,y), 



(7a) 
(7b) 
(7c) 
(7d) 



Using the Boltzmann relations [H + ] = p s exp[— Q(x, y, 0)] 
and [B~] = p s exp[+$(x, y, 0)] for the ion densities at 
the surface z = 0, we can express the associative and 
dissociative surface charge densities a± = [SH + ], [SB - ] 
and al = [S d 



], respectively, as 



m A (x,y) 



l + K A /p s exp(±$(a;,y,0))' 

m%(x,y) 

l+p./^exp(±$(x,y,0))' 



(8a) 



(8b) 



for both types of charging. Here, $(x, y, 0) is the surface 
potential, which for heterogeneous surfaces is a function 
of the lateral coordinates x and y. Note that all surface 
chemistry is encoded in m± D {x,y) and equilibrium 

constants K± D , which we consider as input parameters 
in this work. 

For simplicity we only consider stripe-like surface 
inhomogeneities in this work, such that the charge and 
the surface potential only depend on x and not on y, 
see Fig. [T] Translational invariance in the y-direction 

follows directly by assuming m±/ D (x,y) = m± (x). 
Note that the potential $(x, y, z) — $>(x, z) does depend 
on a lateral coordinate and a normal one. 



The common case that only a small fraction of the 
sites charges (<r -C m) necessarily corresponds to the 
case that the unity in the denominators of Eqs. 



and (8b I is negligible compared to the other term. In 



other words, for all cases that the surface is only charged 
to a fraction of its maximum it is safe to omit the unity 
from Eqs. (8a) and (8b). For most associative reactions 



this approximation is fine, as was already argued by Gra- 
hame in 1947. [3T] For dissociating surfaces one should be 
more cautious, especially at low ionic strengths p s . 32] In 
the remainder of this study we employ this assumption 
throughout, as it reduces the number of independent pa- 
rameters vastly. This is immediately seen by introducing 
the dimensionless total surface charge density 



y(x) = AttXbk- 1 (a A (x) + a°(x) - a A (x) - a°{x)) , 



which reduces Eq. ( 8a ) and ( 8b ) within the small-charge- 
fraction limit to 

y(x) = F + (x)cxp(-$(x,0)) - Y-(x)exp(+$(x,0)), 

(10) 

where 



Y ± (x) 



[(x) 



} {x)K. 



2K A 



(11) 



Note that Eq. ( 11 ) defines a dimensionless quantity which 



is not necessarily of the order of unity. Depending on the 
material properties and ion concentrations, both Y + {x) 
and Y~{x) can easily vary over many decades. From 
now on we refer to Y + {x) and Y~(x) as the positive and 
negative chargeability, respectively. Eqs. ^ and (10) 
give rise to the boundary condition for our model, 



$'(x,0+) = — &(x,Q-) +Y-(x)exp(+$(x,0)) 



y + (x)exp(-$(x,0)). 



(12) 



The latter equation together with Eqs. ([T]) and |3])-([5]) 
forms a closed set of equations to solve the electrostatic 
potential above and within the solid medium. 



B. The homogeneous limit 

To the best of our knowledge, there are no analytic so- 
lutions to the nonlinear Poisson-Boltzmann (PB) equa- 
tion for inhomogeneously charged surfaces. However, it 
can be solved analytically in the case of a single homo- 
geneously charged plate, see for example Ref. [33 . Con- 
sidering a homogeneous surface charge density <j(x, y) = 
<th, the potential can be written as 



( , 1 + 7exp[— kz] 
®h(z) = \ l-7exp[-Kz] 



(13) 



k 1 a H 



where 7 = (y4 + — 2)/y H , and y H = 4n\ 
is the dimensionless surface charge density as defined in 
Eq. @, now for a homogeneously charged plate. Note 
that the solid medium is free of electric fields in this case, 
and because of this, Eq. ( 13 1 is independent of the di- 



electric constant of the plate, ep. By using Eq. (10) for 



(9) 



the mechanism to account for associative and dissociative 
charging we find an expression for the surface potential 
of homogeneously chargeable plates. In general, for a 
(positive) homogeneous chargeability Y + (x) = Yh (and 
Y~(x) = 0) we obtain from Eqs. (10) and (13) that the 
surface potential is given by 
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73 C ° S (l (V^t-- 1 )) if ^ YH < 373' 

$ H (0) = 2hW ?1/3 f9y H +y372nf^) 1/3 .„ , (14) 



[9Y„ + V3yf27Y, 



if Y„ > 



18V3 " '■ti - 373 



The regime Y# < 1 is associated with the regime where 
linear Poisson-Boltzmann theory holds, as we can easily 
check that surface potentials do not exceed unity here. 
By taking the inverse of Eq. ( 13 ) we obtain a relation 



which is well-known from Gouy-Chapman theory. 



yjj = 2 sinh 



*i?(0) 



(15) 



and this can be used in combination with Eq. ( 14 ) to find 



the explicit relation Yh(Yh) between charge and charge- 
ability of a homogeneous plate. The limiting cases are 




if Y H « 1; 
Yl /3 if Y H » 1. 



(16) 



Physically this means that the charging of a weakly 
chargeable surface is relatively efficient, whilst a highly 
chargeable surface only gains charge with the cube root of 
the density of chargeable sites. Both regimes of Eq. ( 16 ) 



can easily be distinguished in a double- logarithmic plot; 
the solid line in Fig. [7] consists of two straight lines with 
slopes 1 and 1/3, respectively, with a cross-over between 
both regimes at Yh ~ 1. 



III. MODE EXPANSION 



that <&(x, z) = <$>(x + L, z) for all x, z, it can be written 



as 



(18) 



where 4> k (z) is the mode amplitude of the Fourier com- 
ponent Tj k (x). One easily checks that functions of the 
form (f) k (z)r]k{x) = c k exp( l -^^)r] k (x), with c k a con- 
stant, solve the PB equation for z < 0, and satisfy bound- 
ary condition d5]). This solution can be put into Eq. pj), 
to yield boundary conditions in terms of the mode am- 
plitudes given by 



e^'(0+) = -^e\ B a k 
lim (f>k(z) = 0. 



27r|fel 



ep^fc(O); 



(19) 
(20) 



The task is now to find PB-like differential equations for 
every mode (j) k (z)r] k (x) in the regime z > by inserting 
the mode expansion ( p~8[ ) into the PB equation ([I]). The 
left hand side can be treated easily, yielding 



V 2 $(r) = ^ U k "{z) 



( 



27rfcY 



M*)]Vk(x), (21) 



We develop a numerical scheme to solve the PB equa- 
tion in the case of a striped patchiness on the surface, 
described by a(x,y) = cr(x), such that the charge dis- 
tribution on the surface only changes in the x, direction. 
We assume inhomogeneities to be periodic with period 
L, such that a{x + L) = cr(x), and we can thus write 



a(x) 



fc= 



M 

^2 cr k r] k (x), 

M 



(17) 



where r] k (x) = —kj- exp( l2 ^ fca ) with k s Z are con- 
veniently normalised (Fourier-) modes and a k are the 
corresponding amplitudes in the expansion, given by 

°fc = K Jq cr (2 ; ) 7 7fc(~ x) Ax. Note that this charge den- 
sity is not known beforehand, as it depends on the sur- 
face potential through Eq. (10). The number M signi- 



fies a high-wavenumber cut-off, that we will empirically 
choose to be large enough to describe the essential large- 
wavelength physics. We will develop a method to calcu- 
late the electrostatic potential $>(x, z) at and above the 
plate, z > 0. Since this is a function which will be sub- 
ject to the L-periodic symmetry as a function of x, such 



where a prime denotes a derivative with respect to z. The 
sinh$(r) on the right side of the PB equation is a non- 
linear function, which gives rise to the couplings between 
mode amplitudes cf> k (z) for all k € Z.[34] This compli- 
cates the calculation of the solution of these modes; in 
general it will not be possible to find solutions for every 
mode separately. 



A. Linear PB 

Only in the case that we describe a weakly charged 
system, such that the potential is small, and thus 
sinh$(r) ps $(r), the PB equation becomes linear, 
V 2 <I>(r) s» « 2 $(r). In this case the modes decouple, and 
using Eq. (21 ) the solution for each mode is found to be 



4>k(z) — a k exp(-K fc z) + b k exp(K fc z), (z > 0) (22) 

where a k and b k are integration constants and K k — 

\J k 2 + (^f^) 2 is the mode-dependent screening parame- 
ter. Note that in the present case of a single plate the 
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coefficient b k vanish because of Eq. (20), and we find by 



applying Eq. ( 19 1 the amplitudes 



Ane\ B cr k L 
eK k L + 27r|fc|ep 



(23) 



Fig. [2^l shows the relation between the screening pa- 
rameters Kk and the (dimensionless) periodicity of the 
system kL. Besides the independence of Ko on kL, it 
shows that K k 3> kq for small kL. This means that 
wave-like inhomogeneities in the electrostatic surface 
potential always vanish within a few wavelengths normal 
to the surface. Therefore, inhomogeneities with short 
wavelengths are screened over shorter distances than a 
Debye length. From Fig. [2^ and from the definition of 
Kk it can be concluded that all inhomogeneities in the 
potential must have essentially disappeared at distances 
of the order of L from the surface as K k ^S> 1/L for all 
|fc|^0. 

Fig. [2}d shows some of the corresponding mode coef- 
ficients afc, for fc = 1,3, and 7, all for ep = (non- 
penetrating fields) and e p = e (index-matched solvent 
and plate). In the figure we divided by clq to show the 
inhomogeneous amplitudes relative to the homogeneous 
background. All coefficients can be calculated by using 
the tTfe following from a surface with charged (y=l) and 
uncharged (y=0) stripes of equal width. This is a fixed 
charge density and therefore we do not account for as- 
sociation/dissociation reactions at the surface here. The 
coefficients a k , and therefore the inhomogeneities in the 
potential, are relatively small but may become large for 
systems with a larger periodicity (kL 3> 1). This is anal- 
ogous to the fact that at small periodicities (with respect 
to k~ x ) the potential is not able to laterally adapt to 
the oscillations in the surface charge and becomes more 
homogeneous, whilst at very high kL it is able to take 
the form of a step function. The choice of the ratio of 
the dielectric constants between the plate and the liq- 
uid is of importance. We see from Fig. [2Jd that if we 
choose an index-matched plate and solvent, correspond- 
ing to the situation that the fields are able to penetrate 
into the solid medium, the coefficients a k become sig- 
nificantly lower (in absolute value) than if this ratio is 
chosen close to zero, where we do not find any fields in 
the medium. In the ep = case, inhomogeneities in the 
charge distribution will have a larger effect on the in- 
homogeneity of the associated electrostatic potential, as 
the polarizability of the solid medium is not able to com- 
pensate any inhomogeneities in the electrostatic potential 
from within the plate. 




Figure 2: The mode-dependent screening parameter n k (a) 
and mode amplitudes a k (b) for several k as a function of 
the periodicity parameter in the system, in (b) for charged 
(y = 1) and uncharged (y = 0) stripes of equal width. The 
solid lines in (b) correspond to the case ep = 0, whilst the 
dashed lines show data for ep — e. 



(exp[$] — exp[— $]) /2 to write the mode expansion of the 
electrostatic potential as 



^ M M 

sinh$(r) = - J] E+(x,z)-- ]J Ej{x,z) 



(24) 



j=-M 



j=-M 



where M is the high-wavenumber cut-off and Ej [x, z) — 
exp [±ej)j(z)r)j(x)] contains the nonlinear dependence of 
the PB equation on (j)j( z )- Using the Taylor expansion 
of the exponential function, and applying that rjj(x) r 



(kL) 



(l-r0/2 



Vn 



j(x), we write Ef(x,z) as a mode sum, 



where high-frequency modes rj n .j for which \n ■ j\ > M 
have been neglected. Note that because of this trunca- 
tion we only need to expand Ej (x, z) up to order M/j 
(and for j — we choose to stop at order M .) As an il- 
lustrative example, E^ (x, z) can be expanded to a mode 
sum as 



B. Nonlinear PB 

In the general nonlinear case, the mode ampli- 
tudes will not be exponential functions of the dis- 
tance like in Eq. ( 22 ) . Instead we now use sinh $ = 



Ef(x,z) = 1 ± (f>i(z)rii(x) 



4>i{z) 2 ri 2 {x) Mz) 3 V3(x) 
6kL 



2VkL 



(25) 



Using the mode-expansion representation of Ej~(x, z), it 
is now a straightforward task to calculate both the prod- 
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ucts of all Ej(x,z) and all E- (x,z), as appearing in 
Eq. (24), and rewrite these products as new mode sums 



with mode amplitudes Vfe(z), 
1 M 

smh$(r) = 2 E (V k + (z)-V k (z))r lk (x) 
k=-M 

+ 0(r) M+ i(x)). 



Note that V±(z) = Vtifai*), M z )> ■ ■ ■ 
that combining Eqs. (26) and (21) yields 



(26) 
(f) M (z)), such 



„2 



v; + ({^W})-y^-({^W}), 



(27) 



for each k, \k\ < M. We solve Eq. (27) for all k it 



eratively as follows. For a given mode-amplitude <p k {z) 
we project out the dependence on all other modes <fij(z), 

j ^ k, by_expanding V^{z) = T,fLo U kA z ) ( l ) k( z ) and 



write Eq. (27) as 



M 



2 M 

yE^W^W)' z>0 - ( 28 ) 



i=0 



As an example we calculate the explicit expressions for 
the monopole (k = 0), 



(±%) 



V kL — 



'-1H>1 
\l kL 



4(kL) 3 /2 



(29) 

where we omitted higher modes |fc| > 1 and left out the 
z dependence for brief notation. If we insert this expres- 



sions into Eq. (28) again, we find 



4>o"vo 



kL 



</> 2 -i0j 
4(kL) 2 



sinh [0o 770] 



which reduces to the planar PB-equation in the case 
that 0±i vanish. Eq. ((28]) together with the BC's (fl9]) 



and ( 20 1 can be used to find the numerical solution to ev- 



ery mode, given the approximations to the solutions for 
other modes. In order to obtain the solution of the full 
problem, we therefore apply an iterative scheme in which 
we go through multiple cycles of solving each PB mode 
equation consecutively, until we find the converged solu- 
tion. This is done by solving for every (f) k (z) the mode 
equations (28) on a 1-dimensional grid. As a starting 
point we use a vanishing solution for all modes. It is 
possible to reduce the computational effort significantly 
by considering systems for which Q(x, z) — $(— x, z), 
such that all 4> k (z) take values in the real space K and 
(j) k {z) = <j)^ k {z). The number of iterations required for 
convergence depends on the degree of nonlinearity in the 
system and is typically in the order of 10 per mode. 



IV. RESULTS 

A. Charging of stripes 

The theory and numerical method we introduced can 
be used to calculate the charging of any periodic config- 
uration of chargeable parallel stripes on a planar surface. 
Here, we focus on a plate coated with stripes with peri- 
odicity L = D 1 + D 2 , where Z?i and D 2 are the widths 
of the striped regions with dimensionless chargeabilities 
Y^ 1 and Y 2 respectively, such that 



Y ± {x) = 



Yi if x £ region 1; 
Y 2 ± if x € region 2. 



(31) 



Regions "1" and "2" are also shown in Fig. [T] We focus 
mostly on = Y, Y{~ = Y 2 = 0, and e P /e < 1, 
which corresponds to uncharged areas (regions "2") 
separating stripes with positive chargeability Y on a 
low-epsilon plate in contact with a high-epsilon liquid 
such as water. For this particular calculation it suffices 
to choose M = 64 modes, a z-grid of N = 2000 points 
extending to z = 5k _1 . 



For a periodicity of 10 Debye lengths, kL = 10, 
a charged stripe of width Di = 3.1k _1 (such 
that D 2 — 6.9k" 1 ), and several chargeabilities 
Y 6 {0.1, 1, 10, 100}, Fig. |3ja) shows the x-dependence 
of the chargeability, which is a step function. This is the 
starting point of our calculation, and we calculate the 
electrostatic potential $(x, z) via the iterative scheme 
that was described above. Fig. [3jb) shows the resulting 
surface potential $(2;, z — 0) for the four values of Y 
of Fig. [3^ a). For each Y, the thin horizontal dotted 
lines show the surface potential in the case that the 
stripe would have infinite width (D\ 3> k -1 ) by using 
the analytical solution to the planar-PB Eqs. ( 14 )-( 15 ), 
i.e. the surface potential of a homogeneously charged 



plate. In the present system the homogeneous surface 
potential limits the actual surface potential from above. 



(30) Since the stripes are relatively broad [D\ > 



), the 

calculated surface potential approaches this limit, which 
we will call the homogeneous limit from now on, at the 
center of the stripe. Typical length scales over which the 
potential varies laterally are clearly of the order K~ , 
as expected. Fig. [3^c) shows the dimensionless surface 
charge y(x), which is related to the surface potential via 
Eq. ( 10 ). We see that the charge density at the center of 



the stripe is well described by that of a homogeneously 
charged plate. By contrast, a relatively high charge 
builds up close to the edges of the stripe. For Y ^> 1, 
in the nonlinear screening regime, the charge density 
at these edges largely exceeds the values we find at the 
center. The reason is the nearby neutral surface, which 
results in a surface potential which is lower at the edge 
of the stripe than at its center. The charging of the 
surface groups, which is normally limited by the induced 
rise of the electrostatic potential according to Eq. (10), 
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Figure 3: Positional (lateral) dependence of (a) the charge- 
ability, (b) the surface potential and (c) the dimensionless 
surface charge density of a striped patchy surface with pe- 
riodicity L = 10k -1 and patch stripe width D\ = 3.1k" 1 , 
for several stripe chargeabilities. The shaded area shows the 
position of the stripe on the plate, and the horizontal dotted 
lines in (b) and (c) show values we would obtain for infinitely 
wide stripes. 



Figure 4: Positional (lateral) dependence of the dimension- 
less surface charge density of a striped patchy surface with 
periodicity L = 10k -1 and stripe width D\ = 3.1ft -1 . We 
indicated the position of the stripe by the shaded areas. Both 
(a) and (b) show a dashed line which shows data correspond- 
ing to Y = 10 as in Fig. [3|c), whilst the solid lines show the 
effect of setting the chargeability of the surrounding plate to 



1 in combination with Y 



combination with Y„ 



0.) 



and (b) Y+ = 50 (in 



can therefore be stronger close to the edges. As a result, 
the average charge density of a stripe can be much 
higher than what one would expect for homogeneously 
charged plates with the same chemical properties. The 
reason that we only observe this effect (deeply) in the 
nonlinear regime is because the surface potentials must 
be significant ($ (x) > 1), such that the Boltzmann 
factors which govern the surface charge distribution, see 
Eq. ( |10p , deviate strongly from unity. 

We checked that the results we obtain do not depend 
on the finite grid size and the number of included 
modes characterized by N and M, respectively; even 
the curves obtained with M — 16 and N = 1000 are 
indistinguishable from all those in Fig. [3j It should be 
noticed that in the case we choose M too small, instead 
of giving inaccurate results, the iterative scheme does 
often not converge anymore, such that no solution is 
found at all. 

The two panels in Fig. [4] each show the charge density 
y(x) for Y = 10 as was calculated in Fig. |3jc) for 



kDi =3.1 and kL — 10 by a dashed line, while the full 
curve denotes y(x) in the case of a modified parameter 
set. It illustrates cases where the surrounding surface 
is chargeable as well. Fig. Qa) shows a situation where 
the surrounding stripe is able to charge up slightly 
negatively (with chargeability — Y). The presence of 
strong peaks of (opposite) excess charge at the edges of 
the stripes demonstrates that this presence of a charge- 
able surface with opposite sign of charge enhances the 
charging of the adjoining area. If the surrounding stripe 
is chosen to be more positively chargeable (F 2 + = 50) 
than the original one, as in Fig. |4]^b), the stripe and 
the surrounding surface change roles. At the interface 
between the two regions the potential now is higher than 
at the center of the stripe, having a value somewhere 
in between the homogeneous limits of both stripes. 
This causes the charge density to peak just outside the 
original stripe whilst inside it shows sharp minima at 
the edges. 

The distribution of charge at the chargeable stripe 
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Figure 5: Surface charge density profiles y{x) of a striped 
patchy surface with Y = 10 for various stripe widths Di and 
spatial periodicities L, here as a function of the distance to 
the center of the charged stripe. In view of the symmetry we 
only plot half a period. The data in (a) is calculated using 
the fixed value of kL = 10, whilst in (b) we use L = 2D\ . The 
vertical dashed lines indicate the edges of the stripes, beyond 
which y(x) = 0. 



depends on the width and the spatial periodicity of the 
stripes. This is depicted in Fig. [5j where the charge 
density y(x) for Y = 10 is plotted like in Fig. |3jc), now 
for various D\ and L. Figs. (5^a) and (b) show the effect 
of changing the width of the stripes but deal differently 
with the size of region in between the stripes. In Fig. [5^ a) 
the stripe periodicity is kept fixed at L = 10« -1 , such 
that the stripe-fraction of the surface increases with 
increasing stripe width D\. In Fig. [BJb) the stripe 
periodicity is set to L = 2D\, for several D\, such that 
the surface coverage fraction of the stripes remains fixed 
at 50%. The numerical results in Fig. [5^ a) clearly show 
that smaller stripes at a fixed stripe periodicity L gain a 
higher charge density. From Eq. ( 10 ) it follows that the 



maximum charge density will occur for infinitesimally 
thin stripes, since the surface potential at the stripe will 
vanish in that case. In this limit, the corresponding 
charge density is y(x) = Y, for all x at the stripe. 
The results for a constant surface coverage in Fig. [5jb) 
fraction also show that thin stripes gain the highest 
average charge density, although the maximal charge 
density at the edges decreases a little. The increase is 
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Figure 6: Results of the numerical calculation for the elec- 
trostatic potential of a striped patchy surface with period- 
icity L — 5k _1 , stripe width D\ — 2n~ l and chargeability 
Y = 100, showing the absolute value of the first 4 mode am- 
plitudes of the electrostatic potential in (a), and a contour 
plot of the electrostatic potential at both sides of the inter- 
face using the calculated mode amplitudes (b). A contour 
line borders the darkest area at $(x,z) > 2.75 and the light- 
est area at &(x, z) < 0.5, with steps of 0.25 in between. The 
vertical dashed line in both figures denotes the location of 
the interface between the solid medium (left) and the solvent 
(right). For these calculations we choose an index-matching 
plate and liquid, e — ep, such that the electric field is able to 
penetrate into the solid medium. 



less pronounced than in Fig. [5]ja) since now the width of 
the charge-neutral region in between the stripes scales 
with the width of the stripes. Therefore, thin stripes in 
Fig. [5^b) are relatively close to their neighboring stripes, 
and will hinder each other in gaining charge. One 
can show that that the maximum charge density will 
again be found for infinitesimally thin stripes, for which 
the system reduces to an essentially homogeneously 
chargeable plate with chargeability Yh — YD\/L. 

Fig. [6] shows the electrostatic potential around the 
plate for kL = 5, kD x = 2, = Y~ = 0, Y+ = 100 
and ep = e. This choice of the dielectric constants is 
such that electric fields do not vanish inside the solid 
medium. In Fig. [6]ja) the logarithmic plot shows nearly 
exponential decay of the mode amplitudes <fik(z) for 
k = 0, . . . , 3. Deviations from a straight line are due to 
nonlinear couplings between the modes, and slopes far 
away from the surface are the mode-dependent screening 
parameters k^- It can be seen that, for this choice of 
ep, the nonzero modes also give rise to an electric field 
inside the solid medium, and for high k the screening 
on both sides of the interface is equally efficient, as 
there are almost no ions involved. Fig. [6^b) shows a 
contour plot of the potential $(x, z) around the charged 
interface. It shows a local increase of the potential close 
to the stripe, which is caused by a local high charge 
density. The inhomogeneity in the potential persists 
for a few screening lengths into the liquid, as the stripe 
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relatively much space in between. This is intuitively 
clear since for stripes that are rather separated the 
charging of the edges is not hindered by the charge 
of neighboring stripes. Moreover, narrow stripes have 
relatively more edge surface. For broad stripes this 
gain at the edges is small compared to the charge at 
the center, and we see from the line for kDi =3.1 and 
kL = 10 that the charge density indeed approaches the 
homogeneous limit (Z?i = L), given in Eqs. ( 14 1-( 15 ) 
and limiting behaviour ( |16[ ). By contrast, the narrower 
stripe of width Di = 0.1k -1 and the same periodicity 
L = 10k -1 has, for Y 3> 1, a charge density that is a 
factor ~ 2 higher than the homogeneous limit, while the 
relative increase of y compared to yn can be a factor 
5 for extremely narrow patches with D\ — 0. 0031k -1 
and L = 0.1k -1 , as can be seen in the inset of Fig. [7] 
Furthermore, we note that y = Y in the low-charge limit 
Y <C 1. This is not surprising either, as in this limit 
any charge, regardless its position on the surface, is 
essentially "alone" in a vanishingly small potential; the 
vicinity of a neutral area "2" or an essentially neutral 
stripe "1" is indistinguishable in that case. 



Figure 7: Average charge density on a chargeable stripe sur- 
rounded by charge-neutral surface for various stripe widths 
D\ and periodicity L, denoted in units of K~ , plotted as a 
function of the stripe chargeability Y. The medium has a 
dielectric constant ep < e. The dashed line corresponds to 
the infinite stripe-width limit for which we use the analyt- 
ical result for homogeneously chargeable plates. Using the 
same data, the inset shows the charge density relative to this 
homogeneous limit. 



width is larger than the screening length, and thus 
ionic screening is the dominant type of screening there. 
In the solid medium ionic screening is not possible, 
and the inhomogeneity persists over a range of the or- 
der of the width of the stripe as was found below Eq. ( 18 1 . 



We now return to the case of chargeable stripes on an 
otherwise neutral plate, and define the average charge 
density on a stripe 



1 

A. Jx=0 



y(x) dx. 



(32) 



Fig. [7] shows y as a function of chargeability Y for 
multiple choices of the stripe geometry parameters 
D\ and £>2- All curves show, as expected, that the 
charge increases with the chargeability. It illustrates, 
however, that the charging of the stripes strongly 
depends on their width and the distance in between two 
successive stripes. Charge densities are increased cither 
by narrower stripes or by larger stripe-stripe distances, 
the strongest charging occurs for narrow stripes with 



B. Analytical approximation 

As a way to better understand and quantify the 
charging mechanism of the patches, we now propose a 
method to estimate the stripe charge densities by an an- 
alytic procedure and check if the right physics emerges. 
The key element in the method is the fact that the 
charging of different surface groups is correlated. This 
correlation is mediated by the electrostatic potential 
and logically the (longest) correlation length is of the 
order of the screening length, as long as we use the 
assumption that ep <C e. This correlation implies that 
a small area of charge-neutral plate around every stripe 
is involved in the charging as well. For the geometry 
of current interest this means we can think of an extra 
strip of width s on each side of the charged patch of 
width Di, as illustrated in Fig. [HJ Now we presume 
that the stripe including the extra area charges up like 
a homogeneously charged plate, such that we can apply 
the analytic expressions ( 14 )-( 15 1 from planar PB theory 
with an effective (decreased) chargeability Y — Y — — — 
for the enlarged stripe. 



D ± +2s 



The strip width s depends on the local screening 
length, and therefore we set s = min(aK -1 , D2/2), 
with a a fit parameter and R~ l the effective screening 
length defined below. The minimum condition is used 
to prevent overlaps, as the size of the additional strip 
area cannot exceed the size of the uncharged region in 
between the charged stripes. The effective screening 
parameter R is determined by the average electrostatic 
potential around the edge. The relation stems from 
PB-linearization procedures and is given for example 
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Figure 8: Pictorial representation of the method used to 
estimate the charging of narrow stripes. The darker areas 
adjoining the shaded stripes resemble strips of charge neutral 
surface assumed to be involved in the charging as well. 



in ref. [35H37] . Since the local surface potential at the 
strip must be somewhere in between zero and the value 
for a homogeneously charged plate with chargeability 
Y, we estimate the effective screeni ng length to be their 
average, such that k" 1 = k -1 / y/cosh (0 + (0))/2. 
The procedure is now to use this effective screening 
length in the calculation of the estimated strip width s to 
obtain the effective chargeability Y. In the homogeneous 
limit this gives a charge density y for the stripes plus 
the side strips via Eqs. (151 and (13|.The average charge 



density y on the original stripe follows by assigning the 
charge density to a smaller surface, y 



y_Dl+2S 



For the fit parameter a — 3/8, Fig. |9]shows the result- 
ing y as a function of Y for several kL and kDi , together 
with the numerically determined "exact" results based on 
Eq. ( 32 1 . The agreement between the two is good, within 



a maximal error of 25%, and we find the right trend 
with Y which suggests that this approximation holds for 
even higher stripe chargeabilities as well. It is expected 
that the fits will remain good or even become better for 
larger stripe distances or stripe widths than investigated 
here, since for broader stripes the relative amount of edge 
surface decreases. Moreover, the choice of parameters 
kL = 10 in combination with kD\ < 3.6 describes a 
system in which the stripes are separated by multiple 
screening lengths, which is close to the kDi oo limit 
where the patches do not mutually interact anymore. 



V. CONCLUSION 

Although the importance of charge inhomogeneities 
on the interactions between surfaces was already men- 
tioned in many other studies, the relation between 
chemical inhomogeneity on the surface and the resulting 
inhomogeneity in surface charge was, to our knowledge, 
never studied in detail. We developed a numerical 
method, based on expansion into Fourier modes, to find 
the electrostatic potential and charge distributions close 



100 




Figure 9: Average charge density on a chargeable stripe sur- 
rounded by charge-neutral surface, for various stripe widths 
D\ and separations D2, plotted as a function of the stripe 
chargeability Y^. The medium bares a dielectric constant 
tp < t. The solid lines show results from the approximative 
method as described in the text, whilst the dotted lines show 
results from the numerical calculations. 



to an interface with an inhomogeneous distribution of 
chargeable chemical surface groups. We focussed on 
the case where the chemical groups are clustered in 
stripe-like regions, with different chargeabilities, either 
of the same or different sign. For the case of striped 
chargeable patches on an otherwise neutral surface we 
showed a significant increase of charging of surface 
groups as the width of the stripes decreases below the 
Debye length. For very small patches this increase of 
the surface charge density can easily be an order of 
magnitude, although patches that have little spacing 
show significantly less increase of charging because 
neighboring patches hinder each other's charging. From 
the numerical results we arrive at the observation 
that the edges of the patches are able to charge up 
optimally and will contribute significantly to the total 
charge of small patches. Our calculated numerical 
results for the average charge density on a patch are 
in fine agreement with a very basic analytical model 
(and one fit parameter) which employs an effective 
patch size in combination with results for homogeneous 
plates. Our results are stepping stone towards the 
study of interactions between heterogeneously charged 
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surfaces, where charged patches may induce charges onto 
nearby (otherwise) neutral surfaces, thereby generating 
nontrivial forces and torques. This could be relevant 
for understanding self-assembly of patchy nanoparti- 
cles or proteins. Work along these lines is in progress. [3"5] 
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